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ABSTRACT 

We argue that the velocity dispersions and masses of galactic bulges and spheroids are 
byproducts of the feedback that regulates rapid black hole growth in protogalaxies. 
We suggest that the feedback energy liberated by accretion must pass through the 
accreting material, in an energy-conserving flux close-in and a momentum-conserving 
flux further out. If the inflowing gas dominates the gravitational potential outside the 
Bondi radius, feedback from Eddington-limited accretion drives the density profile of 
the gas to that of a singular isothermal sphere. We find that the velocity dispersion 
associated with the isothermal potential, a, increases with time as the black hole 
mass M grows, in such a way that M cx a 4 . The coefficient of this proportionality 
depends on the radius at which the flow switches from energy conserving to momentum 
conserving, and gives the observed M — a relation if the transition occurs at ~ 100 
Schwarzschild radii. We associate this transition with radiative cooling and show that 
bremsstrahlung, strongly boosted by inverse Compton scattering in a two-temperature 
(T p ^> T e ) plasma, leads to a transition at the desired radius. 

According to this picture, bulge masses Mb are insensitive to the virial masses of 
their dark matter haloes, but correlate linearly with black hole mass. Our analytic 
model also explains the M& — a (Faber-Jackson) relation as a relic of black hole 
accretion. The model naturally explains why the M — a relation has less scatter than 
either the M — Mb (Magorrian) or the Faber-Jackson relation. It suggests that the 
M — a relation could extend down to very low velocity dispersions, and predicts that 
the relation should not evolve with redshift. 

Key words: galaxies: formation — galaxies: nuclei — galaxies: bulges — accretion, 
accretion disks — quasars: general — black hole physics 



1 INTRODUCTION 

Models offering explanations for the "M — cr" relation — the 
tight correlation between the mass M of the central black 
hole in a galactic nucleus and the velocity dispersion a of 
the galaxy's bulge or spheroid (Ferrarese & Merritt 2000; 
Gebhardt at al. 2000; Tremaine et al. 2002) — fall into two 
broad categories. Global models typically appeal to feed- 
back associated with the growth of the black hole. Kinetic 
energy (Silk & Rees 1998; Blandford 1999), radiation en- 
ergy (Sazonov et al. 2005), bulk momentum (Fabian 1999; 
King 2003), or radiation pressure (Fabian, Wilman & Craw- 
ford 2002; Murray, Quataert & Thompson 2005) liberated 
by accretion reacts back on the infalling material over large 
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distances, stopping both black hole accretion and star for- 
mation in the galaxy when the black hole reaches a certain 
size. Focusing on the response of the galaxy on scales well 
outside the black hole's gravitational sphere of influence, 
these models have paid relatively little attention to the de- 
tails of the accretion process itself. In contrast, local models 
relate the growth of the black hole to details of its immedi- 
ate nuclear environment (e.g., Zhao, Haehnelt & Rees 2002; 
Adams et al. 2003; Miralda-Escude & Kollmeier 2005), but 
typically have focused on dynamical processes rather than 
feedback effects. 

Recently di Matteo, Springel & Hernquist (2005) have 
linked local and global feedback effects, using computer sim- 
ulations of merging galaxies (see also Springel, di Matteo 
& Hernquist 2005). Although their phenomenological re- 
alizations of star formation, radiative cooling in a multi- 
phase medium, and black hole accretion and feedback are 
extremely crude, they are able to reproduce the M — a rela- 
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tion remarkably well. Interestingly, they find that the M — a 
correlation is insensitive to the gas fraction in their model 
galaxies, and that the details of star formation and super- 
nova feedback also have little effect on the black hole mass. 
These results suggest that the feedback regulating black hole 
growth operates on local scales — near the black hole's ac- 
cretion radius or closer in — rather than solely on the global 
scales usually considered. The insensitivity to gas fraction 
presumably results because the gas mass is somehow "maxi- 
mized" on the scales where the accretion rate is determined; 
likewise, on these scales black hole feedback is far more im- 
portant than that due to stars. 

In this paper, we present a simple analytic model for the 
local effects of feedback due to growth of a black hole in a 
galactic nucleus. When the gas supply is plentiful, the black 
hole is able to accrete at (or slightly above) the Edding- 
ton limit, but feedback modifies the gas flow pattern out to 
large distances, far beyond the Bondi radius tb = GM/a 2 . 
If the gas dominates the gravitational potential at r > tb, 
feedback forces the density distribution to approximate a 
singular isothermal sphere, p cc r~ 2 , with a unique veloc- 
ity dispersion that satisfies a oc M 1 ^ 4 . We show that the 
constant of proportionality, which depends on fundamental 
physical constants and processes occurring close to the black 
hole, is consistent with the observed M — a relation under 
plausible assumptions about the gas microphysics. Remark- 
ably, we find that the M — a relation is satisfied not just for 
the final black hole mass, but at every stage during black 
hole growth: a increases along with M. 

The plan of the paper is as follows. In §2 we discuss 
the assumptions of the model and derive the main features 
of the gas flow pattern. Because of the simple nature of the 
model — we invoke turbulent and transsonic energy trans- 
port processes using simple scaling laws — we do not at- 
tempt to derive very detailed results. We find that the flow 
depends on a transition radius between energy-conserving 
feedback, at small radii, and momentum-conserving feed- 
back at large radii. We associate this transition with radia- 
tive cooling, and in § 3 we demonstrate that the flow is likely 
to cool via strongly Comptonized bremsstrahlung at ~ 100 
Schwarzschild radii, under plausible assumptions about the 
gas microphysics. The fact that much of the feedback energy 
goes into radiation, albeit at large radii, limits to feedback 
energy flux to roughly the Eddington limit. In §4 we discuss 
the implications of the model for rapid black hole growth at 
close to the Eddington limit, arguing that a should evolve 
with M in such a way as to preserve M oc a 4 at all times. 
We consider the growth of the bulge or spheroidal compo- 
nent of the galaxy, as a relic of feedback-regulated accretion, 
in §5. We argue that the bulge mass (as opposed to the ve- 
locity dispersion) is regulated by the potential of the dark 
matter halo. However, the properties of the bulge depend 
very weakly on the virial mass of the halo, and the resulting 
bulge mass Mb turns out to scale roughly linearly with the 
black hole mass: M b ~ 10 3 M, consistent with observations 
(Magorrian et al. 1998). The M — a and M - M b relations, 
together, reproduce the Faber-Jackson (1976) relation link- 
ing M b and a. We summarize our main results in §6, where 
we also discuss several predictions and further implications 
for galaxy formation theory. 



2 FEEDBACK MODEL 

Suppose an accreting black hole of mass M liberates energy 
at a rate L f = fMc 2 = £L E (M), where L E (M) is the Ed- 
dington limit corresponding to mass M. We suppose that the 
accretion flow is radiatively inefficient close to the black hole; 
therefore, one might expect the flow to resemble an ADIOS 
(Adiabatic Inflow-Outflow Solution: Blandford & Begelman 
1999, 2004) and to develop a powerful outflow. However, we 
also assume that the incoming gas has sufficiently low an- 
gular momentum that it can be regarded as roughly spher- 
ically symmetric. We posit that under such conditions the 
surrounding gas stifles the wind, i.e., the liberated mechan- 
ical energy does not establish a separate channel through 
which to escape to large distances before interacting with 
the ambient gas. Instead, the power that would normally 
drive the ADIOS wind must pass through the inflowing gas 
on its way out. 

We envisage two ways for the accreting gas to carry 
the feedback energy flux without being blown away. The 
first can be regarded as a form of convection. Let Lf ~ 
4-Kpvfr 2 , where v t is a characteristic eddy speed (assumed 
to be the same as the transport speed) for the energy flux. 
The characteristic energy transport speed cannot exceed the 
Keplerian speed, vk = c(r/r 9 ) -1 ^ 2 , without blowing apart 
the inflow. Adopting this estimate for v t , we find that the 
density obeys 

p(r) ~ x i , (i) 

TgOT 

where m p is the mass of a proton, ot is the Thomson 
cross section, r g = GM/c 2 is the gravitational radius of 
the black hole and x = r/r g . This is analogous to the re- 
sult obtained by Gruzinov (2002) , for non-radiative spherical 
accretion with phenomenological (transsonic) heat conduc- 
tion. It is also analogous to the density law in the CDAF 
(Convection-Dominated Accretion Flow) model for rotat- 
ing, radiatively inefficient flows dominated by ion pressure 
(Narayan, Igumenshchev & Abramowicz 2000; Quataert & 
Gruzinov 2000). In both cases the accreting gas is required to 
carry a conserved outward energy flux by mechanical means. 

Defining the net inflow speed by 4,-Kpv- lu r 2 — M, we have 
fin ~ f~ 1 cx~ 3 /' 2 . For / not too small, v ln quickly becomes 
much smaller than vt- The inflow is very subsonic, nearly 
stopped by the energy flux it is forced to carry. From the con- 
dition of hydrostatic equilibrium, the gas pressure must obey 
p ~ pGM/r oc r~ 3/2 . We then find that p/p 5/3 oc r~ 2/3 , i.e., 
the specific entropy is a strongly decreasing function of r. 
This means that the inflow is strongly convectively unstable. 
The outward heat flux can be regarded as a very inefficient, 
highly saturated convective flux, making it reasonable to as- 
sume that v t is close to the maximum possible speed, vk- 

The second mode of energy transport is relevant where 
radiative cooling is important. Under the conditions we en- 
visage this occurs in the outer parts of the flow, outside a 
"cooling radius" r c . At r > r c much of the convective en- 
ergy flux is radiated away, but the associated momentum 
flux, Hf(r c ) ~ Lf/v t (r c ), is retained since radiation car- 
ries away relatively little momentum. Because of the rapid 
cooling, gas pressure alone cannot support the flow against 
free-fall collapse. To maintain near-hydrostatic equilibrium, 
the momentum flux must therefore satisfy 
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~ -4n (l - JL) p(r) 5 (r)r 2 , (2) 

where <;(r) is the local gravitational acceleration and L ~ 
£Lb is the radiative luminosity liberated at r > r c . We vi- 
sualize 11/ (r) as an internal flux of momentum, carried by 
predominently radial turbulent motions, rather than as an 
outflowing wind. 

The factor 1 — L/Le at in eq. (|5J implies that the feed- 
back energy flux cannot exceed the Eddington limit, even 
though the flow is assumed to be radiatively inefficient close 
to the black hole. Since we are interested in rapid accre- 
tion in the gas-rich environment of a protogalaxy, we will 
assume that the flow is Eddington-limited and that I « 1. 
Since Hf(r c ) ~ 4irp(r c )g(r c )r 3 , according to eq. © the mo- 
mentum flux should be approximately conserved at r > r c . 

We define the Bondi radius r_g as the place where the 
gravitational acceleration due to the black hole equals g ga ,i(r) 
due to the stars, gas, and dark matter in the galaxy (but not 
including the black hole). At r > Tb the Eddington limit 
associated with the enclosed mass exceeds Le{M), hence 
we can set 1 — L/Le « 1 in eq. P). Integrating eq. @ from 
rs to some larger radius r, we obtain 

II/ (r) ~II/(r c ) -4tt j pgr 2 dr, (3) 

J r B 

with g ~ pgai(r). We identify three qualitatively distinct 
cases. If pgr 3 increases as a power of r, the momentum flux 
will be quickly depleted and the feedback will fail. Mat- 
ter will slump toward the center until pgr 3 is roughly inde- 
pendent of r. On the other hand, if p decreases so rapidly 
that dllf/dr ~ 0, there will be enough momentum flux to 
drive a strong outflow at r > rs- Once the local mass sup- 
ply is depleted the accretion rate, and therefore the feed- 
back, will decline, allowing matter to fall toward the black 
hole and increase p. We assume that the third case, inter- 
mediate between the first two, is the one that is chosen 
by nature. We conjecture that p(r) adjusts itself so that 
p{r)g{r)r 3 ~ Hf(r c ) ~ constant (to within a factor ~ lnr). 
This self-regulation is the crux of our feedback argument; 
its validity will have to be checked by more detailed means 
— presumably through numerical simulations. In § 4 we will 
explore its ramifications, but first we discuss the value of the 
cooling radius, r c . 



3 COOLING RADIUS 

As we shall see in §4, the coefficient of the predicted 
M — a relation depends on the radius at which the feedback 
makes a transition from energy- conserving to momentum- 
conserving, and that our predictions match the observa- 
tions for x c ~ 100. In this section we show that such val- 
ues of x c are plausible, provided that two conditions are 
met: 1) the electrons and ions are coupled thermally by 
Coulomb scattering and are permitted to have different tem- 
peratures; and 2) the dominant cooling mechanism is ther- 
mal bremsstrahlung, strongly amplified by inverse Compton 
scattering. 

We adopt the density profile given by eq. (Q with 1=1 
and assume that the protons track the virial temperature, 
kT p /m p c 2 ~ x^ 1 . Defining a dimensionless electron tem- 
perature 8 e = kT e /meC 2 and dividing volume heating and 



cooling rates by m p c 3 / (r 2 ar) , we obtain the dimensionless 
Coulomb energy transfer rate from protons to electrons, 

c pc ~ o.oi3 ( x- 1 - —e e ) e: 3 '^- 1 (4) 
V m p J 

(Stepney & Guilbert 1983, where we have assumed a 
Coulomb logarithm A ~ 20). Our justification for assuming 
strong Comptonization comes from the fact that the flow 
at r < r c is very optically thick to electron scattering, with 
a characteristic optical depth t(x) ~ x 1 ^ 2 , and Coulomb 
coupling is likely strong enough to maintain an electron 
temperature 8 e > 0.1. Thus, the Compton "y— parameter" , 
y w 46> e (l + 46» e )r 2 (Rybicki & Lightman 1979), is large 
and Comptonization is highly saturated, with low-frequency 
photons boosted into a Wien peak at the inverse Compton 
temperature (which equals the electron temperature, under 
these conditions) . The dimensionless bremsstrahlung cooling 
rate is given by 

Ccbr ~ 2.9 x 10~ 8 .4 el /2 x~\ (5) 

where A is the amplification factor due to Comptoniza- 
tion of the bremsstrahlung, discussed at length by Ry- 
bicki & Lightman (1978: see their eq. [7.74] and preced- 
ing discussion). The cooling radius is defined by the con- 
dition that the volume-integrated bremsstrahlung emissiv- 
ity equal Lj; this is approximately equivalent to the lo- 
cal condition C c br(x c ) = 3x^ 3 - Simultaneously demanding 
Cpc = C c br, we can solve for x c and the electron temperature 
8 e (x c ). The result is weakly dependent on M — applying 
the prescription in Rybicki & Lightman (1978) we obtain: 
x c = 74, 68, 63, 58; 6 e {x c ) = 0.46, 0.44, 0.42, 0.40 for 
M = 10 4 , 10 6 , 10 s , 10 10 M Q , respectively. The Compton 
amplification factor ranges from 275 to 487, over the same 
range of masses. 

At r < r c , the electrons are nearly isothermal at 9 e (x c ), 
while at larger radii the continued importance of Comp- 
tonization and bremsstrahlung emission (which creates a 
large number of photons at low energies, even if the to- 
tal emissivity is declining) is likely to lead to a decreasing 
Compton temperature for some range of r > r c . The optical 
depths and densities in the flow are so large that expulsion 
of the gas at large distances by X-ray pre-heating is unlikely 
(Sazonov, Ostriker & Sunyaev 2004; Sazonov et al. 2005). 
For example, the ionization parameter £ = L/nr 2 is < 10 at 
r > r B - 

Cooling by Comptonized thermal cyclotron emission is 
unlikely to be as important as bremsstrahlung under the as- 
sumptions of our model. To estimate the contribution from 
cyclotron losses, we assume that the magnetic energy density 
is ~ 10% of the ion pressure. The emission at low harmon- 
ics is heavily self-absorbed, and we find that the dominant 
source of seed photons for Comptonization comes from emis- 
sion at harmonics n cyc < 75 (Bekefi 1966; Takahara & Tsu- 
ruta 1982). The Comptonized emissivity scales with radius 
roughly cx x~ 3 , implying that cyclotron emission either cools 
the flow close to the black hole, or not at all. Seed photon 
emission from high harmonics increases the Comptonized 
emissivity by a factor ~ n 2 yc ; nevertheless, we find that un- 
der optimistic conditions, cooling by Comptonized cyclotron 
emission fails by a factor of a few. We note that synchrotron 
radiation and Compton scattering by highly relativistic elec- 
trons may also contribute to radiative losses, but since they 
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deplete only the relativistic (nonthermal) electron popula- 
tion they are unlikely to be decisive in cooling the flow. 

Therefore, we conclude that it is plausible that 
Eddington-limited flows dominated by mechanical feedback 
can remain radiatively inefficient out to x ~ 60. In the re- 
mainder of this paper we will treat x c as a parameter with 
a fiducial value of 100, allowing readers to gauge the sensi- 
tivity of our predictions to details of the gas microphysics. 



4 RAPID BLACK HOLE GROWTH AND THE 
M-a RELATION 

Now let us consider the rapid growth of a black hole during 
the epoch of galaxy formation. We assume that there is so 
much gas falling into the central regions of the protogalaxy 
that the gas mass dominates locally over that of both the 
dark matter and the stars that have formed so far. At r > rs 
we have 

ffgai(r) ~ i-nGpr, (6) 

where the factor 4-zr anticipates the radial dependence of p. 
The only thing that can limit the infall of gas is the feed- 
back from the black hole, but at the same time increasing 
the gas density at tb will increase the black hole accretion 
rate. As noted in § 2, we assume that the total feedback en- 
ergy flux reaches a plateau at roughly the Eddington limit. 
This limiting feedback applies even if the accretion rate is 
much larger than the nominal limiting value of ~ IQLe/c 2 
(where we assume an accretion efficiency of ~ 0.1): in this 
case, most of the liberated binding energy would be advected 
into the black hole (Begelman 1979) and / would be -C 1. 
Substituting for p(r)g ga i(r)r 3 ~ 11/ (r c ) from above, we ob- 
tain 

Thus, feedback due to Eddington-limited growth of the black 
hole causes the inflowing gas at r > rs to distribute itself 
as a singular isothermal sphere. Note also that p does not 
depend on the thermal state of the accreting medium, i.e., 
on the details of heating and cooling processes at r > rs. 

The velocity dispersion associated with the gas distri- 
bution is given by a 2 ~ 2ixGpr 2 tx M 1/2 . Thus, M oc cr 4 at 
all stages during the Eddington-limited growth of the black 
hole. Quantitatively, we can write the proportionality in the 
forms 




~2.4xl0 8 (^) 4M 01 (8) 

where 0200 = o /200 km s _1 . We note that this relationship is 
similar but not identical to the expressions derived by King 
(2003) and by Fabian (1999). This similarity should not be 
a surprise, since the feedback effect in all three models relies 
on momentum balance. There are two fundamental differ- 
ences between our model and the earlier ones. First, in our 
model the black hole mass determines the velocity dispersion 
of the bulge, rather than determining the bulge mass. As we 
will see below, the bulge mass is determined by the potential 
of the dark matter halo, once the bulge velocity dispersion 



is fixed by black hole feedback. Second, black hole and bulge 
growth in our model do not stop because the infalling gas is 
somehow "blown away" by the feedback. A given dark mat- 
ter halo could, in principle, host black holes and bulges with 
a wide range of masses and velocity dispersions (although, in 
practice, the supply of gas available to build both bulge and 
black hole is probably correlated with the mass of the halo) . 
What our model implies is that M, a, and (as we shall see) 
the bulge mass are all correlated, whatever their final values. 
Considering the crude approximations used in our analysis, 
the theoretical correlation © is remarkably close to the ob- 
served M — a relation (Tremaine et al. 2002), provided that 
the transition radius r c is not too different from the value 
estimated in §3 based on bremsstrahlung cooling. The sym- 
bolic expression is intended to display the dependence of the 
M — a relation on fundamental parameters. In particular, 
the first quantity in parentheses is roughly the square of the 
ratio of the classical electron radius, ro = e 2 /m e c 2 ~ cr^/ 2 , 
to the gravitational radius of a proton, Grn p /c 2 . 

As the gas supply is depleted, for whatever reason — 
star formation, mass loss from the outer parts of the sys- 
tem, and/or simply the limited gas content of the halo — 
the accretion rate drops below the Eddington limit and the 
black hole mass levels off. The final mass of the black hole 
determines the final value of a associated with the gas distri- 
bution and the stellar spheroid that eventually forms from 
it. Thus, our simple feedback mechanism explains the ob- 
served M — a relation as a relic of the co- evolution of M 
and a. 



5 BULGE MASS AND THE FABER-JACKSON 
RELATION 

Although our feedback model relates a to M, it does not 
directly determine the mass of the bulge that grows in tan- 
dem with the black hole. At radii where gas dominates the 
gravitational potential, the enclosed bulge mass increases as 
Mb{< r) tx pr 3 ~ M(r /rs). Were it not for the existence of 
a dark matter halo, this behavior could extend out to very 
large radii (limited mainly by the availability of gas) and the 
bulge could be many orders of magnitude more massive than 
the hole. However, if the dark matter density distribution is 
shallower than that of the gas, as it is believed to be in the 
central regions of galaxies (Navarro, Frenk & White 1997, 
hereafter NFW), it will eventually dominate the potential. 
Where the inner NFW density distribution, pnfw oc r~ l , 
dominates the potential, we have p oc r~ 3 from the momen- 
tum conservation condition. The enclosed bulge mass thus 
goes from increasing linearly with r to increasing logarith- 
mically. The radius where Pnfw = p thus determines the 
mass of the bulge. 

To evaluate the bulge mass predicted by our model, we 
assume the "concordance" ACDM model with (at z = 0) 
O m = 0.3, D.a = 0.7, and h = 0.65. The inner NFW density 
profile is then given by 

PNPw(r) » 9 x 10- 24 M° ;i ° 2 7 r£ c [£(z)] 2/3 g cm" 3 (9) 

(Komatsu & Seljak 2001), where M vir = 1O 12 M v ,i 2 M is 
the virial mass of the halo, £(z) = [fl m /n m {z)](18TT 2 + 82x- 
39x 2 )/100 with x = fl m (z)-l (Bryan & Norman 1998), and 
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we have used c « l2.8M-° 2 13 h' 013 (l + z)" 1 (Bullock et al. 
2001) for the concentration parameter. 

Equating pnfw to p we obtain the bulge radius, 

n ~ 7.4 (^) V4 M-tf [ 5 (z) r 2/3 kpc, (10) 

where Ms = M/10 8 Mq, and the ratio of bulge mass to black 
hole mass, 

M b r b ( X c X 1 / 2 .,-0.07 < t , -.,-2/3 /„x 

where we have used eq. JHJ. Thus, we find an approximate 
proportionality between the mass of the black hole and the 
mass of the bulge, roughly consistent with the Magorrian et 
al. (1998) relation. The virial mass of the dark matter halo 
enters only to the -0.07 power. The M — a relation then gives 
Mb tx a 4 , recovering the Faber-Jackson (1976) relation for 
spheroids. 

Although the M — a relation is independent of redshift, 
the ratio between M and the bulge mass Mb can exhibit 
weak z— dependence due to the factor [£(z)]" 2 / 3 . For the 
concordance model with a cosmological constant, £(z) de- 
creases from ~ 1 at z = to 187r 2 f2 m /100 w 0.53 as z — > oo. 
For z < 0.5, £(z) « (1 + z)" 0,8 . The general trend is for 
Mb/M to increase slightly with redshift, but by less than a 
factor of 2. For ACDM, M b /M is larger by about 30% at 
z = 0.5, 40% at z = 1, and 50% at z > 5. This modest 
trend will be further weakened by the fact that black holes 
observed at a given z formed over a range of redshifts. 



6 DISCUSSION AND CONCLUSIONS 

We have described a simple feedback model relating the 
Eddington-limited growth of supermassive black holes in 
young galaxies to the properties of the inflowing gas. The 
model requires that the feedback fluxes of energy and mo- 
mentum pass through the material surrounding the black 
hole, on all scales except perhaps in the immediate vicinity 
of the event horizon. This implies that the incoming gas not 
hit a centrifugal barrier until approaching the black hole. 
Although the sub-Keplerian angular momentum of the flow 
probably reflects boundary conditions in the outer parts, we 
speculate that the turbulent fluxes of energy and momentum 
sustain quasi-spherical flow at smaller radii by transporting 
angular momentum outward. 

The model predicts that, once the inflowing gas domi- 
nates the gravitational potential in the galaxy's core, feed- 
back forces the gas to develop the density profile of a singular 
isothermal sphere, with a velocity dispersion that depends 
on the instantaneous mass of the black hole. In contrast to 
most other feedback models, which assume that the black 
hole grows to a certain limiting mass in a fixed galactic po- 
tential, our model implies that the velocity dispersion asso- 
ciated with the bulge, a, increases as M increases, in such a 
way that M oc a 4 at all times. 

Our model reproduces three observed correlations: the 
M — a relation, the Magorrian (M — Mb) relation, and the 
Faber-Jackson (Mb — cr) relation. If it is correct, the model 
implies that spheroids form as the result of self-regulated 
accretion onto supermassive black holes and grow in tandem 
with them. Of the three correlations, the M — a relation is 



the most direct, and therefore it should be no surprise that it 
seems to be the tightest observationally. Both the Magorrian 
and Faber-Jackson correlations also involve the properties 
of the dark matter halo, and therefore should exhibit more 
scatter. 

Our feedback model makes some clear predictions and 
has a number of interesting consequences for galaxy for- 
mation. Since the coefficient of proportionality between M 
and a 4 depends firstly on the maximum feedback luminosity 
(which we assume is close to the Eddington limit) and sec- 
ondly on gas microphysics close to the black hole (where 
electron-ion coupling via Coulomb scattering and Comp- 
tonized bremsstrahlung determine the cooling radius) , there 
should be no limit to the range of velocity dispersions over 
which the M — a relation can apply. The only requirements 
are that the growth of the black hole be dominated by 
Eddington-limited accretion and that the feedback operate 
according to our scheme. (These conditions could be vio- 
lated, for example, if the inflowing gas were rapidly rotating 
at large radii or if mergers of smaller black holes had played 
a large role in black hole growth.) In particular, the M oc a 4 
scaling could extend to extremely low velocity dispersions, 
where cooling rates at tb differ dramatically from that of 
bremsstrahlung. 

The model predicts that there should be no dependence 
of the M — a relation on redshift. However, the ratio between 
the bulge mass Mb and M could show a weak trend with z, 
depending on cosmological parameters. For the concordance 
ACDM model with a cosmological constant the trend should 
be toward higher ratios Mb/M with increasing z. Black holes 
and bulges forming at z = 0.5 could have Mb/M ~ 30% 
higher than black holes and bulges forming at z ~ 0; Mb/M 
could be as much as 50% higher for black holes forming at 
z > 5, but the amplitude of the effect would be significantly 
reduced if the black holes observed at z D b s formed over a 
range of redshifts much larger than z bs- 

Although the ratio M/Mb is well-determined by the 
model, the final mass of the black hole, and of the bulge, 
could be sensitive to processes that trigger the rapid inflow 
of gas into the center of the halo. These may depend on the 
details of galaxy mergers or tidal encounters, and lead to a 
large scatter in, e.g., ratios of black hole or bulge mass to 
halo mass or circular velocity. This could explain why black 
hole masses appear to correlate with bulge properties but 
not with the properties of disks (Kormendy 2001). However, 
it is likely that the supply of gas available for growing the 
bulge and black hole is to some extent correlated with the 
mass of the halo. A linear correlation appears to be neces- 
sary in order to explain the entropy excess of galaxy groups 
and clusters in terms of the energy injected by a central 
black hole (Roychowdhury et al. 2004). 

The model also predicts that rapidly accreting black 
holes should be heavily obscured, with a characteristic col- 
umn density N ~ 1.5 X 10 25 (x c /100) 1/2 cm" 2 , independent 
of M. Black holes accreting at close to the Eddington limit 
are therefore expected to be Compton-thick, as previously 
discussed by Fabian (1999; see also Fabian et al. 2002) and 
Hopkins et al. (2005). 

The feedback mechanism we have described should also 
be relevant under certain conditions when the gas outside tb 
is not self-gravitating and/or the black hole is accreting at 
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less than the Eddington limit. Possible consequences for the 
fueling of active galactic nuclei will be discussed elsewhere. 

Feedback mechanisms like the one we have described, 
in which the feedback energy must pass through the in- 
flowing matter in a (presumably) highly turbulent fashion, 
present a significant challenge to numerical modelers. From 
the perspective of state-of-the-art models like that recently 
described by di Matteo et al. (2005: see also Springel et al. 
2005), our feedback process represents "subgrid" physics. 
However, the basic physical elements of feedback, radiative 
cooling, and conserved momentum and/or energy fluxes are 
present in both treatments, and it is not impossible that 
some of the effects we describe are implicit in their simula- 
tions. 
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